Nano-optomechanics with localized carbon-nanotube excitons 
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We analyze a framework for optical manipulation of the motional state of a suspended carbon 
nanotube based on deformation potential exciton-phonon coupling. The excitons are confined using 
an inhomogeneous axial electric field which generates optically active quantum dots with a level 
spacing in the milli-electronvolt range. A transverse field induces a tunable parametric coupling 
between the quantum dot and the fiexural modes of the nanotube. We show that this interaction 
enables efficient optical ground-state cooling of the fundamental mode and will allow access to 
quantum signatures in its motion. 



Optical transducers underpin a host of high preci- 
sion measurement techniques and recent developments 
in optomechanics suggest that they may enable quan- 
tum limited control of a macroscopic mechanical de- 
gree of freedom pQ. Given the versatileness of mechan- 
ical non-linearities this would provide an alternative to 
atomic systems for fundamental tests of quantum me- 
chanics and the development of quantum technologies 
[2]. Paradigmatic goals in this direction are the prepa- 
ration of a mechanical resonator in its quantum ground 
state [3j [U [5j [6] and the demonstration of quantum sig- 
natures in its dynamics [7] [8] . These endeavors are seri- 
ously hampered by the mechanical quality of typical ma- 
terials. In this respect suspended single-walled carbon 
nanotubes (CNTs) [9] are emerging as a unique candi- 
date. Indeed recent transport experiments in these sys- 
tems have demonstrated strong coupling of charge to vi- 
brational resonances [lOj [11] and ultra-high mechanical 
frequency quality-factor (fQ) products [T2] . These de- 
velopments raise the question of which are the prospects 
for optical manipulation of motional degrees of freedom 
in CNTs. The standard paradigm in optomechanics is 
based on an optical cavity whose frequency is modulated 
by the motion of one of its mirrors or of a dielectric 
object inside it via radiation pressure effects pQ. How- 
ever this approach becomes inefficient for resonators with 
deep subwavelength dimensions and low polarizabilities 
like CNTs. 

Here we propose a solution to this conundrum based on 
an alternative way of inducing coherent optomechanical 
transduction which exploits the unique properties of ex- 
citons in semiconducting CNTs [I3j [I4j [15] [16] . The role 
of the optical cavity is played by an excitonic resonance 
of the CNT that couples parametrically to the motion via 
deformation potential electron-phonon interactions [T7] . 
Homodyne detection of the output field of the two level 
emitter afforded by the excitonic resonance allows then 
to perform a continuous measurement of the mechani- 
cal amplitude. This procedure, which could be imple- 
mented using the differential transmission technique [30] , 
is analogous to cavity- assisted schemes pQ and equiva- 
lent to ion-trap measurements. We envisage a suspended 



CNT where the center of mass (CM) of the exciton is 
localized via the spatial modulation of the Stark-shift 
induced by a static inhomogeneous electric field. We an- 
alyze a tip electrode configuration that effectively engi- 
neers a pair of tunable optically active nanotube quan- 
tum dots (NTQDs) with excitonic level spacing in the 
meV range corresponding to a confinement length below 
lOnm. The quantum confinement is induced by the in- 
homogeneity in the field component along the CNT axis 
En . In turn the normal component E± can be used to in- 
duce a tunable parametric coupling between the exciton 
and the fiexural motion of the CNT. This allows for op- 
tical ground-state cooling of the fundamental mode at an 
ambient temperature in the Kelvin range. A major ad- 
vantage of this mechanical resonator-NTQD system with 
respect to prior scenarios [4] is the possibility of realiz- 
ing a mechanical analogue of the strong-coupling regime 
of cavity-QED [18] with an "optomechanical coupling" 
in the 100MHz range. Furthermore, this coherent cou- 
pling can be switched on and off on demand which offers 
rich possibilities for deploying quantum-optical schemes 
to demonstrate quantum signatures in the motion [19]. 

The electronic structure of a semiconducting CNT can 
be understood in terms of graphene rolled into a cylin- 
der [20]. In the absence of a magnetic field there is a 
single bright level: the singlet bonding direct exciton 
\KK*) + \K' K'*) (the conjugated wavefunctions corre- 
spond to the hole) [21]. Threading a small Aharonov- 
Bohm flux (/>ab renders the antibonding state \KK*) — 
\K f K f *) weakly allowed so that its spontaneous emission 
rate T can be tuned. Thus we focus on the En direct 
excitons \KK*) ± \K'K'*), whose zero-field splitting lies 
in the meV range [15], and consider their deformation 
potential (DP) coupling to the low frequency phonons 
corresponding to the compressional (stretching) and fiex- 
ural (bending) branches for </>ab = [33 . To obtain a 
tractable model for the excitonic wavefunction suitable 
for analyzing exciton CM confinement and the exciton- 
soft phonon coupling we adopt: (i) the k • p graphene 
zone-folded scheme following Ref.[2T]but neglecting inter- 
subband transitions [22], and (ii) an envelope function 
approximation within each subband. For the latter we 
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adopt the parameterization developed in Ref.|22]but take 
the Bloch function at K, K' as determined by (i) and the 
assumption of electron-hole symmetry. This leads to the 
following singlet direct exciton wavefunctions: 

IVw±> =\ (\K n , + K^_) ® \F nm ) ± \K'_ m , + K'\_) 

®\K m ))®(\n) + \im (i) 

where the envelope functions \F nm ), \F' nra ) satisfy 
F nm (z e ,z h ) = Fn m (z h ,z e ) — note that the inter-valley 
mixing preserves the total momentum and angular mo- 
mentum. The indices n, m correspond to the quanti- 
zation of the single particle azimuthal momenta (n = 
0, ±1, ±2, . . .) so that n = m = for En excitons. The 
associated subband electronic ID Bloch functions are 
given by \K ni ±) = [e^-'M* /Vfo](sgn[3n - v\, ±1, 0, 0) 
and |i^ ± )' = [e^ n + l/ / 3 )^/v / 2^](0,0,sgn[3n + v\,±l) 
where + and — label, respectively, the conduction and 
valence bands, v — ±1 corresponds to the "type" of semi- 
conducting CNT [22:, and cp is the azimuthal angular co- 
ordinate (cf. Fig. [I]). Here the unit vectors correspond 
to graphene 2D Bloch functions with wavevector K [K'\: 
(1, 0, 0, 0) [(0, 0, 1, 0)] is constructed from p z orbitals cen- 
tered at the A sites with phase e~ %e l 2 \e %e l 2 \ at the cen- 
tral site, while (0,1,0,0) [(0,0,0,1)] has the p z orbitals 
centered at the B sites with phase e %e l 2 \—e~ %e l 2 \ at the 
central site [20] where is the CNT's chiral angle. 

The symmetric tip electrode configuration sketched in 
Fig. [I] with voltages V\ and V2 allows independent tun- 
ing of E± and E\\ as the reflection symmetries imply 
that they are determined respectively by (V\ — V<i)l2 and 
(Vi +V2V2. For the parameters that allow to confine the 
exciton's center of mass (CM), E±(z) and E\\(z) can be 
regarded as constant across the CNT's cross section and 
the length scale over which they vary appreciably is much 
larger than the excitonic Bohr radius. It follows that for 
sufficiently weak magnitudes (see below): (i) the effect 
of E\\ is dominated by intrasubband virtual transitions 
whose effect on the CM motion can be treated adiabati- 
cally, while (ii) E±_ leads to a weak perturbation ex e ±llf 
that only induces intersubband virtual transitions. 

More precisely, within each pair of subbands n, m (i) 
can be described by H[ nt = ez' ^(zqm) + 0[z' 2 } where 
z' = z e - z h , z C m = [m M z e + m( m h h ]/[m^ + wS m \ 
and E\\(z) = — ^(z). We consider E\\(z) much weaker 
than the critical field to ionize the exciton so that 
(i/; nrri ± I ifi n t I i/j n m±) is much smaller than the binding en- 
ergy. The latter allows to derive an effective Hamilto- 
nian for the ground state manifold of the quasi- ID hy- 
drogenic series associated to n, m by adiabatic elimina- 
tion of the corresponding excited manifolds. This effec- 
tive Hamiltonian for the CM motion has a potential part 
whose leading contribution is second order in i^i nt and 
yields the effective confining potential: V^ m ^(zcM) = 
— a^^E 2 (zqm)- In a classical picture the field polarizes 
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FIG. 1: (color online) Top right: schematic view of the device 
(both the tip electrodes and the CNT are suspended). Top 
left: electric field magnitude and field lines as calculated with 
COMSOL for Vi = -4.5V and V 2 = 5.5V assuming for sim- 
plicity cylindrical symmetry around the tips' axis. Bottom: 
Perpendicular field component and effective confining poten- 
tial along the CNT axis. The resulting exciton CM proba- 
bility distribution is also shown (a.u.). The tailored confine- 
ment generates a QD molecule (double well) with minima at 
\z\ = zo comparable to the distance between the tip electrodes 
[E'll(^o) = 9.6V/i _1 and inter- well tunneling is negligible]. We 
consider the restrictions imposed on the fields by the need 
to avoid dielectric breakdown and sizable static deflections 
(CNT polarizabilities are estimated from Refs. 123]) . 



the exciton that thus experiences a force proportional 
to the gradient of the squared field. The characteristic 
level spacing of the hydrogenic series implies that the 

excitonic polarizability satisfies a^ m ^ ~ 2e n _ m (0)a^ m ^ 

where cr^ m ^ is the exciton Bohr radius and e n - m (q) an 
appropriate dielectric function. In particular eo (0) = ey~ 7 
corresponds to the intrinsic permittivity along the CNT 
[24] . As shown in Fig. [l] we have calculated the field 
generated by this electrode configuration for typical pa- 
rameters using FEM. These result for a (9,4) tube [22] . 
in a zero point motion (Jq^ ~ 9nm [34] which yields a 
level spacing AE = 2.4meV. Given the subwavelength 
separation between the two wells, for each doublet only 
the bonding state is bright so that the degeneracy is ir- 
relevant. Thus the wavefunctions for relevant E\\ are well 
approximated by Eq. ([T]) with modified F nm (z / ^ z cm) = 
f eh^\ z/ )f cm \ z Cm) which can be taken to be Gaussian 
for the ground state (henceforth |^oo±)) [22]. 

We now consider the effect of E± on |^oo±) to lowest 
order in perturbation theory. In principle the linear cor- 
rection |^oo±) involves contributions from all four exci- 
tonic manifolds for which \n\ = 1, m = or n = 0, \m\ = 
1, namely E12, ^21, -E13, -E31. We find that the contri- 
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butions from E12 and E21 vanish identically and obtain 
l4oL) *tEi(\A±,i) + |^°i ± , z ))(F-,o,z|^(4)|Foo) 
where [35] £ = efl/2e ± (En - £13), = -fg(0,z), 

and / labels a complete set of envelope functions for the 
£?i3 and Esi manifolds [36] — e± « 1.6 denotes the in- 
trinsic relative permittivity normal to the CNT axis [24] . 

The Hamiltonian describing the interaction between 
electrons and low frequency phonons has two distinct 
terms: (i) a true DP contribution diagonal in sublattice 
space corresponding to an energy shift of the Dirac point 
and (ii) a bond-length change contribution off-diagonal 
in sublattice space [17]. For the aforementioned exciton 
states we find that the electron-hole and K-K' symme- 
tries imply that finite couplings only arise from (ii): 

H^p =g2 (cos36a x T z — sm36a y ) 

x [Uiptpif) - u zz (f) + 2iu tpz {f)\ . (2) 

Here ?i and &i are, respectively, Pauli matrices in valley 
(K-K') space and in sublattice (A-B) space, and g 2 is the 
off-diagonal DP. Given that the relevant phonon wave- 
lengths A are much larger than the CNT radius R we can 
adopt a continuum shell model [17] with u%j{f) the cor- 
responding Lagrangian strain and use the lowest orders 
in R/X [thin rod elasticity (TRE)]. Both compressional 
and flexural deformations have the structure of a local 
stretching so that the strain components satisfy u ipz ={) 1 

Ucpcp = —cru zz ; with u zz = —R cos (p^if for flexural modes 
and u zz = for compressional modes where (j)f/ c are 
the ID fields* [25] and a « 0.2 is the CNT Poisson ratio 
[26] . Then Eq. ([!]), the aforementioned approximation 
for l^ooL)' and the single particle Hamiltonian ^ allow 
us to obtain the lowest order contributions in the electric 
field to the interaction Hamiltonian i^QD- p h between the 
exciton states |^oo±) and low frequency phonons 

(^ooil^QD-phl^ooi) = -2ug 2 (l + cr) cos3(9 

^2, 



x(F 00 | 



\Foo) (3) 



where we have exploited the completeness of {\F- v qj}}. 

Henceforth, we consider parameters for which the flex- 
ural and compressional branches are expected to present 
resonances with a free spectral range larger than the op- 
tical linewidth of the zero phonon line (ZPL) of the tran- 
sition associated to the state l^oo-)- We focus on laser 
excitation of the latter near resonant with the ensuing 
lowest-frequency flexural phonon red sideband. We con- 
sider a bridge geometry (cf. Fig. [TJ with a length L short 
enough that the relative strength of these phonon side- 
bands is weak [cf. Eq. d5])]. Hereafter, we use the formal- 
ism developed in Ref. [251 and adopt a resonator-bath rep- 
resentation with the resonator mode (annihilation opera- 
tor 60, angular frequency cjo, and quality factor Q) corre- 
sponding to the fundamental in-plane flexural resonance 



that we intend to manipulate and laser cool to the ground 
state. In turn the bath modes include the 3D substrate 
that supports the CNT coupled to the other nanotube 
vibrational resonances. Hence we insert in the effective 
field operators 4>f/ c m Eq. ^ the resonator-bath mode 

decomposition, i.e. ^f(z e ) — > {z e )y/hj 2/jluo (bo + 

&C)) + E d2UX,<1 ( " ^ /On ,.,(n\(h JL h]\ cm A Mc 



dz2 Hze)V h / 2 PMQ)(b q +b\) and ^(z e )- 
Here (j>o{z) is the nor- 



malized resonator ID eigenmode and fi is the linear mass 
density of the CNT, while u XiQ (z) [u z , q (z)] is the x [z] 
component of the CM displacement of the CNT cross 
section at z for the bath mode corresponding to the scat- 
tering eigenmode q and p s is the substrate's density. 

Thus in a polaronic (shifted) representation [4], the 
Hamiltonian for the laser driven NTQD coupled to the 
resonator mode, to the phonon bath (annihilation oper- 
ators b q ), and to the radiation field (annihilation opera- 
tors a,k and couplings g^) reads H = H< 



with (fi = 1): H B = E, 



sys 



q ~q u q u q 

hH.c.) + 



sys T- H{ nt + Hb 

+ E/cK - u L )a\a k , 



9k<r+ 



q 



a z )b q + R.c. (4) 



Here B = e 77 ^ 0-6 ^ S is the laser detuning from the ZPL 
and ft the Rabi frequency, we have introduced Pauli ma- 
trix notation for the optical pseudospin (a z = 1 corre- 
sponds to l^oo-) and a z = —1 to the empty NTQD), 
applied a shift to the phonon modes and adopted a ro- 
tating frame at the laser frequency uol- The parameter 77, 
which characterizes the strength of the exciton-resonator 
coupling (e _r? I 2 is the Frank-Condon factor) is given by 



7] 



~^ 1 + "5|sfey3«»* 



(5) 



where we 



have 



introduced the effective field S± = 
(VL/q^)(F 00 \^(z e )E ± (z e )\F 00 ), q is the TRE phonon 
wavevector for the resonator mode, h = 0.66 A is the 
effective thickness for the continuum shell model [26J, 
£ = lTPa is the CNT Young modulus, and a G = 7.7 x 
10 -7 Kgm -2 is the mass density of graphene. We focus 
on parameters such that 77 < 0.2. Its favorable scaling 
as is a direct consequence of the quadratic flexural 
dispersion. Note that the perturbative treatment of E± 
underpinning Eq. ([5| implies ££_l<C1. Finally, the cou- 
plings ( q and X q to the bath modes lead, respectively, to 
the resonator mode's phonon tunneling dissipation and 
to pure dephasing of the NTQD. The RWA for the ( q is 
justified given their weakness and 77 <C 1. These condi- 
tions and the anharmonicity of the flexural spectrum also 
imply that the flexural X q can be neglected and the pure 
dephasing is dominated by the compressional branch [25J . 

We find that for all environmental couplings the Born- 
Markov approximation is valid and after eliminating the 
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bath phonon modes and the radiation field, we obtain 
a master equation for the NTQD coupled to the res- 
onator with a Hamiltonian contribution given by H sys 
and a dissipative contribution of Lindblad form with col- 
lapse operators: \/f cr_, ^j]j/2a z , ^/uon((Jo)/Qbl ) , and 

y/wo[n((jJo) + 1]/Qb . Here n(uo) is the thermal equi- 
librium occupancy at the ambient temperature and 
is the phonon-induced dephasing rate. Other relevant 
sources of dissipation beyond those considered in Hamil- 
tonian Q can be incorporated by adopting modified val- 
ues of Q \12[ and T [37]. The dephasing rate is de- 
termined by the low frequency behavior of the phonon 
spectral density J(uo) = 7r^ q | | 2 <5(cc; — uo q ) (with q G 
compressional branch). For a bridge geometry the scat- 
tering modes derived in Ref. [25] result in an Ohmic spec- 
tral density J ((J) = 27ra con o;, for uo much smaller than the 
fundamental compressional resonance uo c ^> cjo, which 
naturally leads to = 2na con k B T /h. In turn, it is 
straightforward to determine that the "confined" dimen- 
sionless dissipation parameter satisfies a con = a/i:Q Cl 
where Q c is the clamping-loss limited Q-value of the fun- 
damental compressional resonance [25 and a the dissipa- 
tion parameter that would result for an infinite length. 
The latter can be calculated using Eq. and reads: 
a = 91^(1 + a) 2 cos 2 36/27r 2 hR(Eh) 3 / Y It depends 
on the chirality and may approach unity for small radius 
zigzag tubes. Thus, when the exciton linewidth is dom- 
inated by electron-phonon interactions [27] the phonon 
confinement in our structure will reduce it by at least 
the factor ttQ c ^> 1. Finally, a con <^ 1 warrants the 
Born-Markov approximation in the treatment of the pure 
dephasing in the relevant regime jjj > F/2. 

In complete analogy with the Lamb-Dicke limit we ex- 
pand up to second order the translation operators B and 
adiabatically eliminate the NTQD to obtain a rate equa- 
tion for the populations of the resonator's Fock states. 
This incorporates both the mechanical dissipation and 
the dissipative effects induced by the scattering of laser 
light. The latter result in cooling and heating with rates 
r] 2 A^ that read the same as in Ref. 4 with the quantum 
dot Liouvillian including now the pure dephasing 

7d. As 7£> — > the steady state occupancy for Q — > oo, 
i.e. the quantum backaction limit — be- 

comes independent of ft (in stark contrast to atomic laser 
cooling) and reduces to the same expression valid for the 
cavity-assisted backaction cooling [5] with the cavity de- 
cay rate 1/r replaced by the spontaneous emission rate. 
In the resolved sideband regime and for the optimal de- 
tuning S — — c^o, this fundamental limit yields (r/4c^o) 2 . 

We have studied the temperature dependence of the 
optimal steady state occupancy (6q6 )ss = [ n (^o) + 
r] 2 QA+/uo]/[l + rj 2 Q(A- - A+)/uo] minimized with re- 
spect to S and ft for the following parameters: Q = 2xl0 5 , 
7/ = 0.17, r = cj /3, and h<y D /k B T = 4.7 x 1(T 4 . These val- 
ues correspond, for example (cf. Fig. [I]), to a (9,4) nan- 
otube (a = 0.2) of length L = 120nm {u q /2tt = 725MHz, 



Q c = 908) with £ ± = 147V/i- 1 . For this case an ambient 
temperature T = 4.2K [n(uo) = 120] allows cooling to a 
final occupancy n/ =0.045. In contrast to previous opti- 
cal backaction cooling scenarios [4j [5] there is appreciable 
temperature-dependent pure dephasing that results in a 
steady state occupancy that scales non-linearly with T 
for T>1K. Finally, the corresponding input power for a 
diffraction-limited spot size is below InW rendering ra- 
diation pressure effects completely negligible. 

In the un-shifted representation, for 5 = and after a 
7r/2 rotation of the pseudospin around H sys reduces to 
the Jaynes-Cummings model with the spin degree of free- 
dom afforded by the NTQD states dressed by the laser 
field. The spin-oscillator coupling and resonance condi- 
tion are given, respectively, by g = tjcoq/2 and Q = uoq. 
Thus, given 1/Q <C r) reaching the strong coupling regime 
depends on satisfying t]uoq >T/2. This regime is akin to 
the parametric normal mode splitting in cavity optome- 
chanics [28] and offers a wide range of possibilities for the 
demonstration of quantum signatures in the motion. In 
particular a judicious modulation of r] locked to pulsed 
laser excitation allows to emulate the adiabatic passage 
scheme used in Ref. [19] for performing QND measure- 
ments of the oscillator's energy. This would enable the 
observation of motional quantum jumps. 

In conclusion, we set forth a scheme for optomechanical 
manipulation of nanotube resonators via the deformation 
potential exciton-phonon interaction. This provides a 
high-performance alternative to radiation-pressure based 
schemes [1 for an ultra-low mass and high frequency 
nanoscale resonator leading to large backaction-cooling 
factors and opening a direct route to the quantum be- 
havior of a "macroscopic" mechanical degree of freedom 
0. Most importantly, these breakthroughs rely on a 
lifetime-limited zero phonon line much narrower than the 
smallest CNT linewidths reported so far [13] . Indeed, the 
envisaged NTQDs will allow to suppress the two most 
likely linewidth-broadening mechanisms, namely: inho- 
mogeneous broadening and phonon-induced dephasing 
[27] . by providing a controlled electrostatic environment 
and strong confinement of low frequency phonons. Fur- 
thermore, a doped version of these NTQDs will enable a 
tunable spin-photon interface [29] . 
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